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Abstract 

Matveev and Larkin calculated the parity effect on the ground state energy of 
a small superconducting particle in the regimes where the mean level spacing 
5 is either large or small compared to the bulk gap A. We perform a numerical 
calculation which extends their results to intermediate values of 5/ A. 

I. INTRODUCTION 

Black, Ralph, and TinkhamBi performed tunneling experiments on nanometer scale par- 
ticles of aluminum at temperatures T m 50 mK well below the superconducting transition 
temperature T c ~ 1 K. In some particlesEJ they found a superconducting energy gap larger 
than the level spacing. This gap was slightly enhanced from the bulk gap. In smaller 
particles!, the level spacing exceeded the bulk gap and no evidence of superconductivity was 
found. One may then ask whether all evidence of superconductivity is extinguished in these 
smaller grains. 

Matveev and Larkinll addressed this question within a BCS model 

H = E£j4 a c ja -g £ 4\ c ti°n c i'V (I) 

h a |ej|,|ej/|<w d v ' 

Letting E^ be the ground state energy for the system with m electrons, they defined the 
parity effect parameter 

A P = Ef n+ V - \ (Ef™) + Ef n + 2 )) . (2) 

This parameter measures the degree to which the odd parity ground states have a higher 
energy than the even ones. Matveev and Larkin studied the limit g j 0, to a — > oo, at fixed 
mean level spacing 5 and at fixed value of the bulk gap 

A = 2oo d e- 6/9 . (3) 

They found the asymptotic results 

Ap/A ~ 1 - 5/ A « 1 (4a) 
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A minimum was therefore expected in Ap/A for 8/ A near 1. 

We perform a numerical calculation for Ap/A with constant level spacing 8, and locate 
the position and value of the minimum. We also verify the asymptotic limits of Ap/A are 
given by fl4a[) and ([4b|) . 



II. METHOD 



We consider a collection of states j G J = {—n, . . . ,n} about the Fermi surface and 
study the effective Hamiltonian for these states 

H = 5 E jc+ c ja - g E c/t c ii c i'-l c i'T- (5) 
The effective coupling g is given in perturbation theory by 
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where is the set of integers k satisfying n < \k\ < tOd/8 and /x is the chemical potential. 
We use the effective coupling g in the limit g J. 0, ojd — >• cxd taken at fixed A and fixed 5. 
The coupling ^ can then be written in the form 

9 = Tl (2n+l)5\ ^ 

ln(a n y A ; ) 

where the constant a ni(U ~ 1 for n > 2 and < 5. We nevertheless use the exact a njfl in this 
calculation. 

We are interested in the ground state energy of H for both odd and even electron number. 
In the even case the ground state has no singly-occupied level and can therefore be mapped 
onto the ground state of the spin Hamiltonian 

H = SEi°S-9lL<r}<r? (8) 
for j G J, where the spin operators are defined by 

°t = 44 ( 9 ) 

° 3 = c ii c ;t 

The odd case is similar except for one singly-occupied level (j = 0) at the Fermi surface. 
This level is inert and can be ignored for our purposes. We therefore use the effective 
Hamiltonian (H) with j G J — {0} in the odd case. 

We thus define the Hamiltonians H^ 2n \ H ( - 2n+1 \ and if( 2n+2 ) by (H) and the following 
constraints: 




E<x| = -1 ; j G J 
E<x| = ; j G J-{0} 

E<?| = i ; j e J 



^(2*0 = -5/2 

/i( 2 ™ +1 ) = (10) 



2 



The couplings entering the if( m ) are g even for m = 2n, 2n + 2 and g dd for m — 2n + 1. These 
are given by expression (|7|). As n increases, the ground state energy of H^ m ~) approaches 
that of the original Hamiltonian H with m electrons, up to a constant which is the same for 
all m. 

Since the Hamiltonians if ( 2n ) and if ( 2n+2 ) are related by particle-hole symmetry p, which 
sends o z — > —o z and ^ a T , 

pH^p = H^-g even , (11) 

we need calculate only the ground state energies Ef n ^ and Ef n+l \ We do this by exact 
diagonalization and form the combination 



A D = R(2n+1) — if R(2«) , E (2n+2) 
F 9 2 V 9 ' 9 



E (2n+l) _ E (2n) + 



en 



(12) 
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III. RESULTS 



In Figure 1 we plot the results for intermediate values of S/A. The curves from top 
to bottom are the results for n = 3,4,5, and 6, respectively, and the straight line is the 
strong coupling asymptote. We find excellent convergence in the region near the minimum 
by n = 6. For constant level distribution, we find the minimum Ap/A w 0.7 occurs at 
S/A w 0.9. 

In Figure 2 we show the agreement with the weak coupling asymptote, where 
A p (asymptotic) is given by expression (f4b| ). 

Note. As this manuscript was being completed, there appeared a manuscript by Mastel- 
lone et alQ which obtains similar results. 
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FIGURES 
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FIG. 1. The curves from the top to bottom are for n = 3,4, 5, and 6. They approach the strong 
coupling asymptote indicated by the straight line. 



5 




0.3 1 1 1 1 1 

5 10 15 20 

6/ A 



FIG. 2. The curves from top to bottom are for n = 3,4,5, and 6. All approach the weak 
coupling limit A p (asymptotic) given by expression flibj). 
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